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Abstract 

We present femtosecond transient absorption measurements on vr-conjugated supramolecular 
assemblies in a liigh pump fluence regime. 01igo(;)-phenylenevinylene) monofunctionalized with 
ureido-s-triazine (MOPV) self-assembles into chiral stacks in dodecane solution below 75° C at a 
concentration of 4 x 10^^ M. We observe exciton bimolecular annihilation in MOPV stacks at high 
excitation fluence, indicated by the fluence-dependent decay of l^B„-exciton spectral signatures, 
and by the sub-linear fluence dependence of time- and wavelength-integrated photoluminescence 
(PL) intensity. These two characteristics are much less pronounced in MOPV solution where the 
phase equilibrium is shifted significantly away from supramolecular assembly, slightly below the 
transition temperature. A mesoscopic rate-equation model is applied to extract the bimolecu- 
lar annihilation rate constant from the excitation fluence dependence of transient absorption and 
PL signals. The results demonstrate that the bimolecular annihilation rate is very high with a 
square-root dependence in time. The exciton annihilation results from a combination of fast exci- 
ton diffusion and resonance energy transfer. The supramolecular nanostructures studied here have 
electronic properties that are intermediate between molecular aggregates and polymeric semicon- 
ductors. 

PACS numbers: 78.47.+p, 78.55.-m, 78.55.Kz, 78.67.-n 
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I. INTRODUCTION 



Organic semiconductor materials are now widely used in optoelectronics devices such as 
field-effect transistors,ii*^ light-emitting diodes,i^ and photovoltaic diodes.— Intramolec- 
ular functionality (vr-conjugation) can be tailored with synthetic methodologies , ^1^1 but 
control of intermolecular properties, which are equally important in determining electronic 
properties {e.g. charge transport), is more elusive. Perhaps the most compelling advan- 
tage of many organic semiconductors such as conjugated polymers over molecular and in- 
organic semiconductors is solution processability, making fabrication of elegantly simple 
device structures possible with techniques such as ink-jet printing . ^^1^^ A natural strategy to 
achieve three-dimensional control of intermolecular interactions is to exploit molecular self 
assembly in solution prior to the casting process by means of supramolecular chemistry.^^^ 
With this approach, molecular building blocks are self-assembled to form well-defined, com- 
plex architectures through secondary forces, such as electrophilic/electrophobic interactions. 
These interactions need to be strong enough to lead to spontaneous self-organization but 
weak enough so that the process is reversible. Supramolecular assemblies possess polymeric 
characteristics and are therefore often referred to as "supramolecular polymers". Their 
macroscopic properties can be tuned such that high carrier mobilities, for example, can be 
achieved.— 

Recently, this approach has been applied very successfully to oligo(p-phenylenevinylene) 
derivatives with chiral side-chains and functionalised with ureido-s-triazine, a hydrogen- 
bonding end-group . ^^1^^ It has been shown that these monofunctionalized oligo(p- 
phenylenevinylene) (MOPV) self-organize in chiral stacks in apolar solvents below a tran- 
sition temperature. Here, the secondary interactions used for the supramolecular assembly 
are quadruple hydrogen bonds, tt-tt stacking and solvophobic effects. 

Resonance energy transfer is fundamental to describe exciton dynamics in conjugated 
structures. It has been shown that intermolecular energy transfer is the dominant mechanism 
in polymer films while intramolecular energy transfer determines exciton dynamics in dilute 
polymer solutions . ^^i^^i^^ A quantitative study of intermolecular energy transfer in polymer 
films is complicated by the intrinsic positional and energy disorder of such systems. MOPV, 
on the other hand, is a model system to study energy transfer dynamics as it is possible to 
investigate the effect of intermolecular interactions by comparing optical properties of MOPV 
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in the dissolved phase and the supramolecular assembhes. In supramolecular assembhes, 
the conjugated segments are closely packed in a fashion similar to conjugated segments in 
polymeric films but in a well-defined manner. In a previous paper, we have investigated 
exciton dynamics in MOPV at low excitation fluence with time-resolved photoluminescence 
spectroscopy and we found them to be very similar to those of polymeric films. Fast 
exciton diffusion on chiral supramolecular assemblies was shown to lead to exciton trapping 
and luminescence depolarization. 

It has been documented that at sufficiently high excitation fluence, exciton dynamics 
are dominated by bimolecular annihilation processes. This exciton-exciton interaction has 
been studied in several polymeric system o^'^i^^i^^i^^i^^i^^i^'''i^^i^^i^° (see Gadermaier et alM- for 
a recent review on ultrafast studies of exciton dynamics in polymeric semiconductors). 
Two models have been proposed for the mechanisms of exciton-exciton annihilation and 
have been applied to various polymeric systems. The first model assumes the rate-limiting 
process to be the diffusion and encounter of excitons, with the rate constant being time- 
independent.— >2S*2Si22i2i In the second model, the rate- limiting process is considered to be a 
long-range resonance-energy-transfer between the excitons. The proposed form of the rate 
constant varies but usually has an explicit time dependence due to the process being non- 
Markovian.— i^^jS^iS^iSl These two mechanisms are in competition, and the key parameter 
that determines the importance of the non-Markovian mechanism is the spectral overlap 
of the l^Bu absorption and emission spectra.— The excited-state absorption spectrum is 
generally broad and its position is similar in polyphenylene and polypheny lenevinylene ma- 
terials, but the PL spectral position is more sensitive to chemical composition. As a result, 
exciton bimolecular interactions in green and red-emitting materials tend to be dominated 
by long-range resonance effects, while in blue-emitting materials collisional processes tend 
to dominate bimolecular annihilation dynamics. 

In this paper we describe a comprehensive study of exciton bimolecular annihilation 
in MOPV supramolecular assemblies in dilute solution. Applying femtosecond transient 
absorption spectroscopy, we show that excitation with pulse fiuences above ~ 50/iJ/cm^ 
results in bimolecular annihilation dynamics that dominate on picosecond timescales. This 
is in analogy to solid films of polymeric semiconductors. Long-range exciton bimolecular 
interactions determine the annihilation dynamics as has been reported for related alkoxy- 
substituted polypheny lenevinylene films. These results demonstrate that the electronic prop- 
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Figure 1: Molecular structure of M0PV4. In dodecane, quadrupolar hydrogen bonding of the 
ureido-s-triazine unit leads to formation of dimers, which self-assemble into chiral structures as 
depicted in the schematic diagram. The supramolecular assembly is thermotropically reversible 
with a transition temperature of 75° C at an M0PV4 concentration of 4 x 10~^ M. 

erties of solution-self-assembled supramolecular architectures are similar to those of disor- 
dered polymeric semiconductors. However, supramolecular electronics allow tailored control 
of intermolecular electronic interactions, as well as reversibility of the assembly process. 

II. EXPERIMENTAL 

The synthesis of M0PV4 (see Figure for the molecular structure) has been described 
elsewhere.— The material was dissolved in anhydrous dodecane at a concentration of 4 x 
10~^ M. Measurements were carried in a temperature-controlled absorption cuvette with 
1 mm path length. 

The femtosecond transient absorption apparatus consisted of a laser system based on the 
design of Backus et alM Briefly, 15-fs pulses were produced by a mode-locked Ti:sapphire 
oscillator (KMLabs TS pumped by 4.5 W, 532 nm output from a Spectra-Physics Millennia 
V laser). The pulse train was amplified using a chirped-pulse-amplification scheme at a 
repetition rate of IkHz. The amplifier was a multipass Ti:sapphire amplifier (pumped by 
the 9.5 W, 527 nm output from a Spectra-Physics Evolution X laser). The resulting pulses 



4 



had a pulse width of 100 fs, a pulse energy of 400 /iJ and were centered at a wavelength of 
800nm (1.55eV photon energy). 

The sample was excited with 3.1 eV pump pulses produced by frequency doubling the 
laser fundamental in a 0.5-mm pathlength /5-BBO crystal. The pump power was varied with 
calibrated neutral density filters and the beam was mechanically chopped at 500 Hz. The 
transient absorption of the sample was probed with a single-filament white-light continuum 
(approximately from 1.2 eV to 3.0 eV) produced by focusing a fraction of the fundamental 
into a sapphire window. Probe pulses were delayed with respect to pump pulses with 
a computer-controlled optical delay. All the beams were linearly polarized with parallel 
polarization. The pump and the probe beams were focused on the sample cell to respectively 
~ 50 /im and 100 /xm spots (with an angle between beams of approximatively 20°). The 
probe beam and a reference beam (not transmitted through the sample) were dispersed in 
a 0.25-m spectrometer and detected with a pair of Si photodiodes. Pump-induced variation 
in transmission (AT/T) was extracted using standard lockin techniques. The instrument 
response function had a full width at half maximum of less than 120 fs. 

We applied time-correlated single photon counting (TCSPC) to measure excited-state 
lifetimes. The MOPV solution was excited with a pulsed (20 MHz, 70 ps FWHM, 407 nm) 
diode laser (PicoQuant LDH400). The luminescence was detected with a microchannel plate 
photomultiplier (Hamamatsu) coupled to a spectrometer and TCSPC electronics (Edinburgh 
Instruments Lifespec-ps and VTC900 PCI card). A decay curve was taken at 2.21 eV photon 
energy, and the decay time was extracted from a mono-exponential fit to the data. 

III. RESULTS AND ANALYSIS 

The material investigated here is an oligophenynevinylene consisting of four phenyl rings, 
monofunctionalized with a ureido-s-triazine (M0PV4) to dimerise by hydrogen bonding (see 
Figure P). The absorption spectrum of M0PV4 is shown in Figure |21 together with the pho- 
toluminescence (PL) spectrum at 14°C and 65°C for excitation using femtosecond pulses 
with photon energy just above the peak of the tt-tt* band. A red shift is observed upon cool- 
ing the solution, which has been previously attributed to the formation of supramolecular 
assemblies.~^' ^ Small-angle neutron scattering experiments showed that the average stack 
length is around 160 nm with a radius around 6nm.^^ The stack configuration is not known 
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Figure 2: Absorption spectrum of M0PV4 at 20°C (dotted line), PL spectrum at 14°C (solid line) 
and PL spectrum at 65° C (broken line). The black arrow indicates the excitation photon energy 
in all experiments reported here. 

precisely but according to preliminary quantum-chemical calculations, the oligomer separa- 
tion is approximatively 35 A and the angle of rotation between two neighboring oligomer is 
close to 12°.2i The intermolecular coupling in the stacks modifies the excitonic character- 
istics and leads to the appearance of a red shoulder in the PL spectrum. Above a certain 
transition temperature, which depends on the M0PV4 concentration, the thermodynamic 
equilibrium shifts from supramolecular assembly to a dissolved phase.-^ Within a relatively 
narrow range (± ~ 10°C) about the transition temperature, the positional disorder in the 
stack increases and the average assembly length decreases. At the M0PV4 concentration 
used here, the transition temperature is around ~ 75°C. However, the two PL spectra show 
that at 65°C, the red shoulder has already decreased significantly suggesting that the inter- 
molecular coupling is much smaller at 65°C. Therefore, we can use these two experimental 
conditions (14°C and 65°C) to compare exciton dynamics with and without strong inter- 
molecular coupling. 

Figure El shows the excitation-fiuence dependence of the time- and wavelength-integrated 
PL intensity of M0PV4 at 14°C and 65°C. The dependence is linear with excitation fiuence 
but a sub-linear trend appears at higher fiuence. The transition between the two regimes 
is temperature-dependent. At lower solution temperatures, the PL intensity saturates at 
significantly lower fiuence. Moreover, the deviation from linearity at high excitation fiuence 
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Figure 3: Time- and wavelength- integrated PL intensity at 14° C (filled squares) and at 65° C (open 
circles). The solid and broken lines are the results of fits to the time- and volume-integrated exciton 
population density using the model described in the text. The vertical lines separate linear and 
sub-linear regimes for the data taken at 14° C (solid line) and 65° C (broken line). The two arrows 
point to the excitation fluence of the data used for the stretched exponential fits (see Figure |^. 

is stronger at 14°C. Such behavior is indicative of a time-integrated exciton population which 
develops linearly with low pump fluence and sub-linearly with high fluence. This indicates 
that exciton-exciton annihilation occurs on the supramolecular assemblies at high fluence. 
At low fluence, on the other hand, we are able to probe a regime where exciton-exciton 
annihilation dynamics are negligible. 

Various OPV derivatives have been previously studied using femtosecond transient ab- 
sorption spectroscopy,— i^^iSL^ and three distinct spectral signatures were identified. The 
first is probe-induced stimulated emission (SE) signal overlapped with the PL spectrum 
(around 2.21 eV) while the other responses are photoinduced absorption (PA) signal around 
1.77 eV and 1.46 eV. The SE and the 1.46-eV PA were assigned to l^Bu-exciton dynamics 
while the 1.77-eV PA signal was assigned to polaron absorption. In the following discussion 
we will limit ourselves to the study of the l^B^i-exciton dynamics, as the PA signal due to 
polaron absorption is strongly fluence dependent due to two-step exciton dissociation by 
resonant sequential excitation.— 

Figure |3] shows the time evolution of transient PA signal recorded at a probe photon 
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Figure 4: Absorption transients at a probe photon energy of 1.46 eV and solution temperature of 
14°C (21/uJcm~^ pump fluence, black circles) and 65°C (25/iJcm~^ pump fluence, open squares) 
with fits to the expression (cexp(— — (^)^)- The sohd hne is fit to the data at 14°C (a = 300 ps, 
b = 0.4 and r=1600ps) and the broken hne is the fit to the data at 65°C (a = 300 ps, b = 0.9 and 
r=1600ps). 

energy of 1.46 eV for low excitation fluence (~ 20 /iJcm~^). This corresponds to the regime 
where the PL intensity depends linearly on the excitation fluence (see the arrows in Figure 

and we can therefore neglect exciton-exciton interactions. The transient signal measured 
at 65°C (open squares) decays quasi-exponentially while the decay of the 14°C signal (filled 
circles) displays a much faster initial component. The exponential component is assigned to 
the decay of localized excitons on isolated oligomers. We have shown previously that well 
below the transition temperature, exciton dynamics are dominated by fast diffusion within 
the chiral supramolecular assembly, which leads to trapping and subsequent slower transfer 
of localized excitons.— The ultrafast decay of the signal measured at 14° C can therefore be 
interpreted as a consequence of exciton diffusion along the stacks. 

We have measured time-dependent SE and PA transient signal in M0PV4 solution at 
various pump fiuences, probing at 2.21 eV (Figure El) and 1.46 eV (Figure El), respectively. 
The dynamics of both signals measured at 14°C depend on the excitation fluence, and at 
higher fluence are considerably faster than those measured at 65°C. The fast decay present 
within the flrst 10 ps at 14°C increases in relative contribution and becomes faster as the 




8 




Figure 5: (a) SE signal (probe energy 2.21 eV) at 14°C for 279/iJcm~^ (crosses) and 111 /xJcm^^ 
(circles) pump fluence. The lines are fits with the photophysical model at the corresponding pump 
fluence. (b) SE signal (probe energy 2.21 eV) at 65°C at 836 cm'^ (crosses), 209 ^uJ cm^^ (circles) 
and 70 fiJ cm~^ (squares) pump fluence. The lines are fits to the model at the corresponding pump 
fluence. 

excitation fluence is increased. At those fluences, the initial PA and SE signals (at t = 0) have 
a sub-linear dependence in the excitation fluence (Figure [7j). The deviation from linearity is 
stronger in the data measured at 14°C than for those measured at 65°C. Also, the transition 
between linear and sub-linear regimes at 14°C occurs at lower fluence. These results are 
consistent with the fluence dependence of the integrated PL intensity (Figure Ej). 

The sub-linear excitation fluence dependence displayed in Figures El and [7| are characteris- 
tic of bimolecular exciton annihilation dynamics. Since the transient absorption signal decay 
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Figure 6: (a) PA signals (probe energy 1.46 eV) at 14°C for pump fluence of 557 /iJ cm~^ (crosses), 
209/iJcm~^ (circles), 84/iJcm~^ (squares) and 21^Jcm^^ (triangles), (b) PA signals (probe 
energy 1.46 eV) at 65°C for pump fluence of 557^Jcm^^ (crosses), lll^Jcm^^ (circles) and 
25/xJcm^^ (squares). The lines are the fits to the data at corresponding pump fluence according 
to the model described in the text. 

rates depend on the excitation fluence at 14°C, we correlate the exciton-exciton annihilation 
processes with the formation of supramolecular assemblies. At 14°C, exciton diffusion on 
the supramolecular stacks is fast enough to allow for exciton-exciton interactions at suffi- 
ciently high density. At 65°C, localization of excitons becomes dominant as the disorder is 
increased.— Furthermore, just below the transition temperature, the average stack length is 
shorter as the equilibrium shifts towards the dissolved phased so that the exciton density 
required for exciton-exciton interaction is higher. 
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Figure 7: (a) SE signal (probe energy 2.21 eV) at t = as a function of excitation fluence at 
14°C (black squares) and at 65°C (open circles). This amplitude was determined by fitting the 
transient absorption data in Figures|Hland|niwith a convolution of the instrument response function 
and a multiexponential function, and extrapolating to zero time. This procedure allows accurate 
deconvolution of the instrument response, (b) PA signals (probe energy 1.46 eV) at t = as a 
function of excitation fluence at 14° C (black squares) and at 65° C (open circles), determined as 
described above. The lines through all data are fits to the photophysical model described in the 
text. 

In order to obtain a more quantitative description of exciton bimolecular annihilation dy- 
namics in M0PV4 nanostructures, we have constructed a mesoscopic model that reproduces 
our experimental results and allows us to unravel the mechanism governing the annihilation 
process. The model is based on a rate equation for the exciton population, N, given by 
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The first term is the exciton generation rate, which has the form 



(1) 



G{t) 



aPX {Ngr-S{t))exp{-{t-tof/2a^) 



(2) 



27rr^hc Ngr \jlTxal 



where P is the excitation pulse energy, a is the unsaturated absorbtion coefficient, A is the 
excitation wavelength, r is the Gaussian radius of the excitation beam. Sit) is the number 
of excitons created between to and t (with to the zero-time), A^^^^ is the initial density of 
ground states and Ot is the width of the excitation pulse. To calculate S{t\ we use 



The generation term takes into account the Gaussian distribution in space and time of the 
excitation pulse. The initial density of ground states is taken to be the concentration of 
M0PV4 molecules with the ground state depletion taken into account in each point. 

The second term in equation ^ describes the unimolecular excited-state decay rate, with 
time constant r measured using TCSPC for an M0PV4 solution at 14 and 65°C at a very 
low excitation fiuence. Collisional bimolecular annihilation processes are described by the 
third term, where /? is the diffusion-limited exciton bimolecular annihilation rate constant. 
The general form includes a time-dependence but this term may be neglected if the typical 
diffusion length on experimental time-scales is larger than the extent of the exciton . 
In conjugated systems, this is usually the case; however, this assumption will be discussed 
further in Section HVI 

On the other hand, the long-range time-dependent bimolecular annihilation term, with 
rate constant 7, includes explicit time-dependence. This time-dependence was introduced 
following resonance energy transfer studies and accounts for the dispersive nature of the 
transfer process resulting from the distribution of transfer rates in a random ensemble of 
donors and acceptors, which leads to a non-Markovian depletion of acceptors near the donors. 
The exact time dependence, characterized by the time exponent varies with the spatial 
distribution of acceptors and depends on dimensionality.— However, other factors such as 
disorder may also lead to a power-law dependence of the transfer rate (this point will be 



S{t + dt) = G{t)dt 



(3) 
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discussed further in Section IIV|) and we have therefore considered d in this term as a fitting 
variable. 

The final term in equation Q represents the ultrafast dynamics related to exciton diffusion 
in supramolecular assemblies. We have reported previously dispersive diffusion dynamics of 
excitons in M0PV4 and we will therefore not consider these issues in detail here.«^ To 
take into account effects of diffusion to non-radiative traps, we have added a term in our 
rate equation so that at low excitation fiuence, the solution converges towards the low 
excitation fiuence data. To determine this term, we fit the PA signal measured at the 
lowest excitation fiuence (within the linear intensity regime, see Figure Ej) with a stretched 
exponential taking into account the natural decay (exp(— t/r — (t/a)^)), where a is the 
stretched-exponential time constant and b is the exponent in the stretched exponential. In 
equation ^ e is the amplitude of the contribution of this dispersive process. This method 
follows theoretica l'^^ I'^^i^^ and experiment al^^i^ studies of exciton diffusion to traps were a 
stretched exponential behavior was obtained. The stretched exponential fit to the 65°C PA 
data is very close to an exponential function but the time dependence of the PA data at 14°C 
is close to an exponential of the square root of time. This indicates that diffusion to traps 
and polarization decay influence the PA signal dynamics at 14° C but is less important at 
65°C, as established by our previous studies of exciton dynamics on M0PV4 at low fiuence.— 

The experimental geometry is taken into account in the photophysical model by dividing 
the volume located at the intersection of the probe and the pump beams into small sub- 
volumes with homogeneous excitation fiuence where the rate equation can be solved. The 
measured quantities are then calculated by summing the exciton densities in a mesoscopic 
approach. To calculate the SE and PA signals, the fraction of probe photons transmitted 
through the sample needs to be evaluated. This fraction is related to the exciton density 
and to the SE and PA cross-sections by 

dpi, = ±aNpi,dx (4) 

where p^, is the photon density, a is the cross-section for SE or for PA and x is the sample 
depth. As we are only interested in the sample transmission, the initial probe photon popu- 
lation is arbitrary. The average exciton density at a particular sample depth is determined 
in a similar manner as the total exciton density on the entire volume but the multiplication 
factors needed are the surface fractions excited with each of the particular excitation fiuence 
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Table I: Summary of parameters used in the photophysical model. 





14°C 


65°C 


T fns) 


2940 


2940 


B (cm^ ps~^) 










3.6 xlO"^'^ 





d 


0.5 




e fps~i) 


2.2 xlO~^ 


2.5 xlO"^ 


a, (x)s] 


266 


300 


b 


0.46 


0.9 


a (cm^) 


9.25 xlO"^'^ 


9.25 xlO"^'^ 


Ngr (cm-3) 


2.49 xlO^'^ 


2.49 xlO^'^ 


apA (cm^) 


1.55 xlQ-^^ 


2.5 X 10-16 


asE (cm^) 


4 xlO^i'^ 


1.4 X 10-16 



of the set, at the particular sample depth. 

The parameter r is obtained separately from TCSPC measurements, Ngr from the knowl- 
edge of the M0PV4 concentration, a and therefore a (defined as a /Ngr) from UV-vis ab- 
sorption measurements, and a and b from the stretched-exponential fitting of the low fiuence 
PA signal. To determine the remaining parameters (/3, 7 and d, e, apA and ctse), the PA 
and SE signals at all excitation fiuence, the zero-time signals and the integrated PL data 
are modelled using a global fit. The best results of the fitting procedure are displayed in 
Figures El El inland m with the model parameters summarized in Table HI 

A high bimolecular annihilation rate constant 7 was needed to fit the time-resolved data 
taken at 14°C, but bimolecular annihilation effects were not necessary to reproduce the data 
taken at 65°C. We find that the magnitude of 7 is almost an order of magnitude higher than 
the typical values obtained in polymer films.— i^^i^^*^ Furthermore, the data at 14°C cannot 
be fitted with a time-independent annihilation term, and setting (3 = did not compromise 
the global fits presented here. Variation of the exponent d between 0.4 and 0.7 does not have 
a significant effect on the quality of the fit; however, the best fit is obtained with a square 
root dependence in time [d = 0.5). Figure El shows the exciton density at the maximum 
experimental fiuence used here (~ 840/iJcm-^), as determined by the model. The time- 
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Figure 8: Exciton density at peak excitation fluence against time for the parameters reported in 
TableQlat 14°C (solid line) and 65°C (broken line). Inset: Data plotted over the first 5ps. 

dependent bimolecular annihilation term leads to an ultrafast decay of the exciton density 
at 14°C, but is insignificant at 65°C. Note that the peak exciton population density at both 
temperatures (< 5 x lO^^cm"^) is 1.8% of the ground-state chromophore density used in 
these experiments {Ngr in Table H}, and therefore photophysical saturation is not the origin 
of sublinear fluence dependence observed in Figures El and [7| Such low exciton densities in 
the supramolecular nanostructure point to the need for a high bimolecular annihilation rate 
constant at the photon flux regime used here. 



IV. DISCUSSION 



Our model shows that exciton bimolecular annihilation is necessary to describe excited- 
state dynamics in M0PV4 solution at 14°C while it has a negligible effect at 65°C. To check 
that the high value we obtain for the bimolecular annihilation constant is not a consequence 
of the mesoscopic character of our model, we performed a similar analysis to those devel- 
oped in a range of previous studies.-^i^^^ In this crude analysis, the exciton density was 
modelled with a simplified rate equation with a unimolecular decay term, a time- dependent 
bimolecular annihilation term and a delta function as the generation term. 

_ = Ar„,(()__^,_ (5) 
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exp(-t/r) 



(6) 
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exp(— x' 



For small signals, the PA signal is directly related to the exciton population (AT/T ^ 
—apANx), so we also fitted our data with the analytical solution of the simplified rate 
equation. The exciton population at zero-time, A^o^ was approximated with the known 
density of absorbed photons, which allowed to relate the fitting coefficient to the parameter 
7. The value obtained from the high fiuence PA signals, 7 ^ 2 x 10"^'' cm^ ps"^''^, is very 
similar to the result of our model, which is reproduced by this crude analysis but is closer 
to the M0PV4 physics as it takes into account the variation in excitation fiuence in the 
sample. 

Several studies have interpreted the square-root-of-time dependence of the bimolecular 
annihilation as an effect of the non-Markovian depletion of close acceptor-donor pairs. With 
a point-dipole approximation, the Forster model predicts a time dependence of the excitation 
transfer rate of t~(^"^/^), with D being the dimensionality of the acceptor distribution. This 
result was obtained by applying the method developed by Eisenthal et al. for the three- 
dimensional case to the one- and two-dimensional cases.— If this model is valid to describe 
the exciton-exciton annihilation, the dimensionality characteristic of exciton distribution 
in M0PV4 stacks is three (which corresponds to the best fit, with d = 0.5, but a fractal 
dimensionality between three and two cannot be ruled out). The point-dipole model also 
relates the magnitude of exciton-exciton annihilation rate constant to the Forster radius by 



r is the Gamma function, ^ the exciton lifetime in the absence of transfer and Rq the 
Forster radius. With D = 3, we extract Rq ~ 70 nm with our value of 7, which is un- 
realistic. However, the simple model is based on a random distribution of acceptors in a 
sphere of dimensionality D centered on the donor and it fails to reproduce a distribution of 
excitons in the complex geometry of M0PV4 stacks. We also point out the limitations of 
the applicability of Forster theory in organic semiconductor systems, where the point dipole 
approximation is often unjustified.— 

The coUisional mechanism used in the model does not include any time-dependence. 
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R^ Ti^l^ r (2 - D/Q) 

r(i + D/2) 



(7) 




(8) 
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This is justified if tlie typical diffusion length on experimental time-scales is larger than the 
spatial extent of the exciton. However, the general expression in three dimensions includes 
a time-dependence given by^ 



with D the diffusion coefficient and Rg the collision distance. If we assign the observed 
square-root time dependence to this simple collisional mechanism, we obtain an unrealistic 
value for the diffusion coefficient D ^ 10^ cm^ s~^ with = 1 nm. The collision radius 
Re would have to be around 50 nm for D to be similar to the values found in organic 
systems {D ^ 10^^ cm^ s~^).^^ Furthermore, such high values of D correspond to the time- 
independent regime so that the collisional mechanism cannot explain our data. If we use 
the value of (3 corresponding to the best global fit (which is significantly worse than that 
shown in Figures 01 El El andHj), realistic values for the diffusion coefficient can be obtained 
if -Re > 5 nm.'^'^ The expression used here corresponds to diffusion in a three-dimensional 
lattice, but similar expressions exist for the one and two-dimensions cases. All of these fail 
to reproduce the observed bimolecular exciton annihilation time-dependence and intensity 
dependence. 

With the above considerations, we can conclude that neither the collisional model nor 
the Forster (point dipole) model can uniquely explain the high rate of exciton bimolecular 
annihilation measured in M0PV4 stacks. These two mechanisms represent limiting cases 
where either incoherent exciton hopping or Forster energy transfer are sufficient to describe 
the annihilation. A more realistic description is an intermediate case of fast incoherent 
hopping followed by efficient resonance energy transfer. The square-root-of-time dependence 
of the bimolecular annihilation term extracted from the model in equation is an effect of the 
complex geometry of the supramolecular assemblies and of disorder. Combined with rapid 
exciton diffusion within the stack,^^ bimolecular interactions lead to the high annihilation 
rate. A branching between exciton localization and direct bimolecular annihilation is likely 
at early time. At later times, bimolecular annihilation of localized excitons is likely. This 
picture is analogous to related studies of resonance energy transfer dynamics in blends 
containing small concentrations of M0PV4 in similar oligomers with shorter conjugation 
length.i^ 

A more accurate description of the bimolecular annihilation dynamics would require com- 




(9) 
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paring the time-dependent average exciton density accessed with our measurement and a 
configurationally averaged microscopic description of the transfer processes using a Monte- 
Carlo or a master equation approach. Implementation of these schemes requires a significant 
number of assumptions, including detailed understanding of the structural and geometrical 
disorder. This information is beyond the scope of the present study. 

On the other hand, a detailed description should be based on a quantum-chemical, non- 
Markovian model where the exciton-exciton coupling and electron-phonon coupling are 
explicitly taken into account. OPV systems have been shown to belong to the interme- 
diate coupling regime where the relaxation energy accompanying optical excitation, the 
nearest neighbor exciton-exciton coupling, and the frequency of the most strongly coupled 
intramolecular vibration are comparable.— In M0PV4 supramolecular stacks, the inter- 
molecular exciton-exciton coupling are comparable to typical intramolecular exciton-exciton 
coupling found in polymeric systems.— For such an intermediate case between polymeric 
systems and molecular aggregates that undergo coherent excitation transfer (supertransfer), 
the description of exciton-exciton interaction would need to go beyond a Forster description 
of energy transfer . ^'^i^^ 

We are currently exploring both the microscopic semi-classical and the quantum chemical 
approaches. However, the work presented here demonstrates unambiguously that the exciton 
dynamics in supramolecular nanostructures in solution are similar to those observed in 
polymeric semiconductors, but control of intermolecular electronic interactions is higher 
such that high exciton bimolecular annihilation rates are measured. 

V. CONCLUSION 

We have shown that exciton dynamics in supramolecular assemblies of M0PV4 are very 
similar to the dynamics observed in conjugated polymer films. Supramolecular chemistry is 
a very promising route to construct well-defined semiconductor architectures with polymeric 
properties through controlled reversible assembly. 

We have applied femtosecond transient absorption spectroscopy as a means to explore 
exciton dynamics in a supramolecular nanostructure constructed with a hydrogen-bonded 
oligophenylenevinylene derivative. A simple rate-equation model was applied successfully to 
separate exciton bimolecular annihilation phenomena from unimolecular processes occurring 
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at low fluence.— We have shown that exciton-exciton annihilation does take place on the 
supramolecular assemblies and dominates the exciton djTiamics at high fluence. The mecha- 
nism for bimolecular annihilation is non-Markovian; explicit time dependence is necessary to 
reproduce the experimental data with the model. The rate constant for bimolecular annihi- 
lation is higher than in comparable organic systems and cannot be rationalized with simple 
models used in previous studies. We consider that the electronic properties of this system fall 
in the intermediate electronic coupling regime, where intermolecular electronic coupling is 
comparable to the intramolecular reorganization energy, providing an opportunity to apply 
semiclassical and quantum-chemical methodologies to explore electronic dynamics in this 
promising regime of supramolecular electronics. We are currently pursuing this opportunity. 

Acknowledgments 

The authors thank Prof. Gregory Scholes, Dr. Stefan Meskers, and Dr. David Beljonne 
for fruitful discussions. The work in Cambridge is supported by the UK Engineering and 
Physical Sciences Research Council (EPSRC) and the interdisciplinary research center for 
nanotechnology (Cambridge, UCL, Bristol). CD is an Isaac Newton Scholar. CS is an 
EPSRC Advanced Research Fellow. LMH thanks St John's College, Cambridge for financial 
support. The Work in Eindhoven is supported by the Netherlands Organization for Scientific 
Research (NWO, CW). APHJS is supported with a Fellowship of the Royal Netherlands 
Academy of Arts and Science. The Cambridge-Eindhoven collaboration is supported by the 
European Commission (LAMINATE). 



Corresponding author. Ics271@cam.ac.u^ http : / /www-oe . phy . cam . ac . uk/fast /index . htm| 



^ N. Stutzmann, R. H. Friend, and H. Sirringhaus, Science 299, 1881 (2003). 

2 C. D. Dimitrakopoulos and P. R. L. Malenfant, Adv. Mater. 14, 99 (2002). 

3 R. H. Friend, R. W. Gymer, A. B. Holmes, J. H. Burroughes, R. N. Marks, C. Taliani, D. D. C. 
Bradley, D. A. D. Santos, J. L. Bredas, M. Logdlund, et al.. Nature 397, 121 (1999). 

^ S. Forrest, P. Burrows, and M. Thompson, IEEE Spectrum 37, 29 (2000). 
^ C. J. Brabec, C. Winder, N. S. Sariciftci, J. C. Hummelen, A. Dhanabalan, P. A. van Hal, and 
R. A. J. Janssen, Adv. Funct. Mater. 12, 709 (2002). 

19 



6 B. A. Gregg, J. Phys. Chem. B 107, 4688 (2003). 

^ P. Peumans, A. Yakimov, and S. R. Forrest, J. Appl. Phys. 93, 3693 (2003). 

^ S. Setayesh, A. C. Grimsdale, T. Weil, V. Enkelmann, K. Mullen, F. Meghdadi, E. J. W. List, 

and G. Lcising, J. Am. Chem. Soc. 123, 946 (2001). 
^ C. Ego, D. Marsitzky, S. Becker, J. Y. Zhang, A. C. Grimsdale, K. Miillen, J. D. Mackenzie, 

C. Silva, and R. H. Friend, J. Am. Chem. Soc. 125, 437 (2003). 
^° A. C. Grimsdale, P. Leclere, J. D. Mackenzie, C. Murphy, S. Setayesh, C. Silva, R. H. Friend, 

and K. Mullen, Adv. Funct. Mater. 12, 729 (2002). 

H. Sirringhaus, T. Kawase, R. H. Friend, T. Shimoda, M. Inbasekaran, W. Wu, and E. P. Woo, 
Science 290, 2123 (2000). 
^2 T. Kawase, H. Sirringhaus, R. H. Friend, and T. Shimoda, Adv. Mater. 13, 1601 (2001). 
J.-M. Lehn, Supramolecular Chemistry (VCH, Weinheim, Germany, 1995). 
A. M. van de Craats, J. M. Warman, A. Fechtenkotter, J. D. Brand, M. A. Harbison, and 
K. Mullen, Adv. Mater. 11, 1469 (1999). 

A. P. H. J. Schenning, P. Jonkheijm, E. Peeters, and E. W. Meijer, J. Am. Chem. Soc. 123, 
409 (2001). 

16 E. W. Meijer and A. P. H. J. Schenning, Nature 419, 353 (2002). 

1^ T.-Q. Nguyen, J. Wu, V. Doan, B. J. Schwartz, and S. H. Tolbcrt, Science 288, 652 (2000). 
1^ D. Beljonne, G. Pourtois, C. Silva, E. Hennebicq, L. M. Herz, R. H. Friend, G. D. Scholes, 
S. Setayesh, K. Miillen, and J. L. Bredas, Proc. Nat. Acad. Sci. USA 99, 10982 (2002). 

19 C. Silva, D. M. Russeh, A. S. Dhoot, L. M. Herz, C. Daniel, N. C. Greeham, A. C. Arias, 
S. Setayesh, K. Miillen, and R. H. Friend, J. Phys.: Condens. Matter (2002). 

20 L. M. Herz, C. Daniel, C. Silva, F. J. M. Hocbcn, A. P. H. J. Schenning, E. W. Meijer, R. H. 
Friend, and R. T. Phillips, Phys. Rev. B (accepted). 

21 R. G. Kepler, V. S. Valencia, S. J. Jacobs, and J. J. McNamara, Synth. Met. 78, 227 (1996). 

22 V. I. Klimov, D. W. McBranch, N. N. Barashkov, and J. P. Ferraris, Chem. Phys. Lett. 277, 
109 (1997). 

23 D. Vacar, E. S. Maniloff, D. W. McBranch, and A. J. Heeger, Phys. Rev. B 56, 4573 (1997). 
2"^ E. S. Maniloff, V. I. Klimov, and D. W. McBranch, Phys. Rev. B 56, 1876 (1997). 

25 A. Dogariu, D. Vacar, and A. J. Heeger, Phys. Rev. B 58, 10218 (1998). 

26 G. J. Denton, N. Tessler, M. A. Stevens, and R. H. Friend, Synth. Met. 102, 1008 (1999). 

20 



2^ T. Q. Nguyen, 1. B. Martini, J. Liu, and B. J. Schwartz, J. Phys. Chem. B 104, 237 (2000). 

28 M. A. Stevens, C. Silva, D. M. Russell, and R. H. Friend, Phys. Rev. B 63, 165213 (2001). 

29 M. Shimizu, S. Suto, A. Yamamoto, and T. Goto, Phys. Rev. B 64, 115417 (2001). 

^° C. Silva, A. S. Dhoot, D. M. Russell, M. A. Stevens, A. C. Arias, J. D. MacKcnzie, N. C. 
Greenham, S. Setayesh, K. Mullen, and R. H. Friend, Phys. Rev. B. 64, 125211 (2001). 

C. Gadermaier and G. Lanzani, J. Phys.: Condens. Matter 14, 9785 (2002). 

^2 S. Backus, C. G. D. Ill, M. M. Murnane, and H. C. Kapteyn, Rev. Sci. lustrum. 69, 1207 
(1998). 

P. Jonkhcijm, F. J. M. Hoeben, A. P. H. J. Schenning, and E. W. Meijer (unpublished data). 

D. Beljonne (unpublished data). 

3^ V. I. Klimov, D. W. McBranch, N. Barashkov, and J. Ferraris, Phys. Rev. B 58, 7654 (1998). 
36 D. W. McBranch, J. Kraabel, S. Xu, R. S. Kohlman, V. I. Klimov, D. D. C. Bradley, B. R. 

Hsieh, and M. Rubner, Synthetic Metals 101, 291 (1999). 
3^ B. Kraabel, V. I. Klimov, R. Kohlman, S. Xu, H.-L. Wang, and D. W. McBranch, Phys. Rev. 

B 61, 8501 (2000). 

38 R. C. Powell and Z. G. Soos, J. Lumin. 11, 1 (1975). 

39 T. Forster, Zeitschrift fur Naturforschung 4a, 321 (1949). 

^° K. B. Eiscnthal and S. Sicgel, J. Chcm. Phys. 41, 652 (1964). 

B. Y. Balagurov and V. G. Vaks, Sov. Phys. JETP 38, 968 (1974). 
^2 B. Movaghar, M. Griinewald, B. Ries, H. Bassler, and D. Wiirtz, Phys. Rev. B 33, 5545 (1986). 
^3 B. Mollay, U. Lemmer, R. Kersting, R. F. Mahrt, H. Kurz, H. F. KaufFmann, and H. Bassler, 

Phys. Rev. B 50, 10769 (1994). 

M. Yan, L. J. Rothberg, F. Papadimitrakopoulos, M. E. Galvin, and T. M. Miller, Phys. Rev. 

Lett. 73, 744 (1994). 
^5 L. M. Herz and R. T. Phillips, Phys. Rev. B 61, 13691 (2000). 

G. Cerullo, M. Nisoli, S. Stagira, S. D. Silvestri, G. Lanzani, W. Graupner, E. List, and G. Leis- 

ing, Chem. Phys. Lett. 288, 561 (1998). 
4^ A. Suna, Phys. Rev. B 1, 1716 (1970). 

^8 F. J. M. Hoeben, L. M. Herz, C. Daniel, P. Jonkheijm, A. P. H. J. Schenning, C. Silva, S. C. J. 

Meskers, R. T. Phillips, R. H. Friend, and E. W. Meijer (unpublished data). 
^9 F. C. Spano, J. Chem. Phys. 116, 5877 (2002). 

21 



G. D. Scholes, X. J. Jordanides, and G. R. Fleming, J. Phys. Chem. B 105, 1640 (2001). 
G. D. Scholes, Chem. Phys. 275, 373 (2002). 



22 



